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Abstract 

We study the sensitivity of the ATLAS experiment to Wtb anomalous cou- 
phngs in top pair production with semileptonic decay, pp ^ tt ^ W~^bW~b with 
one of the W bosons decaying leptonically and the other hadronically. Several 
observables are examined, including the W helicity fractions and new quanti- 
ties recently introduced, such as the ratios of helicity fractions and some angular 
asymmetries defined in the W rest frame. The dependence on anomalous cou- 
plings of all these observables has been previously obtained. In this work we 
show that some of the new observables also have smaller systematic uncertainties 
than the helicity fractions, with a dependence on anomalous couplings similar or 
stronger than for helicity fractions. Consequently, their measurement can signif- 
icantly improve the limits on anomalous couplings. Moreover, the most sensitive 
measurements can be combined. In this case, the precision achieved in the deter- 
mination of Wtb anomalous couplings can be of a few percent in the semileptonic 
channel alone. 



1 Introduction 



The three generation structure of the standard model (SM) was completed with the 
discovery of the top quark at Tevatron [1]. Its properties have already been directly 
investigated at colliders [2-5] and, in particular, its mass has been determined to a high 
accuracy [6], better than for any other quark. However, the determination of other 
fundamental properties, like spin and couplings, requires larger top samples, which will 
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be available at LHC. In its first low luminosity phase (10 fb~^/year) LHC will produce 
several millions of top quarks per year and experiment, mainly in pairs through gluon 
fusion gg — > ti and quark-antiquark annihilation qq ti, with a total cross section 
of 833 pb for a top mass rrit = 175 GeV [7,8]. Single top production [9, 10] will also 
occur, dominated by the process bq — > tq' , with an expected cross section of 306 pb [7]. 
Both processes will test the SM predictions for the fundamental properties of the top 
quark, and in particular they will allow us to measure its couplings [9,11-13]. This 
fact is specially important since, from a theoretical point of view, sizeable deviations 
from the SM predictions for top couplings are possible in several SM extensions, as 
for example in supersymmetry [16] and models of dynamical symmetry breaking [17]. 
Indeed, within the SM the Wtb coupling is purely left-handed at the tree level, and 
its size is given by the Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vtt — l- 
But in SM extensions departures from the SM expectation for Vtb are possible [14,15], 
as well as new radiative contributions to the Wtb vertex. These deviations might be 
observed in top pair and single top production at LHC. 

Top pair production takes place mainly through QCD interactions, thus indepen- 
dently of the electroweak Wtb coupling. Additionally, it is likely that the top quark 
almost exclusively decays in the channel t — > W~^b. Therefore, the cross section for 
gg, qq ^ ti W^bW~b is insensitive to the size of the Wtb vertex, as well as to its chi- 
ral and tensorial structure. Still, the angular distributions of top decay products give 
information about the Wtb vertex structure (up to a global multiplicative constant), 
and thus they can be used to probe anomalous top couplings^ In the rest frame of a 
decaying top quark, the energies of the W boson and b quark are fixed by the two-body 
kinematics. Therefore, non-standard Wtb interactions can only influence the following 
groups of observables: 

1. The total width r(t — > Wb), which is very difficult to measure at LHC§| 

2. The helicity fractions of the W boson, which also determine the angular distri- 
butions of its decay products in the W rest frame and their energy distributions 
in the top rest frame. 

^The global normalisation of the Wtb vertex can be determined in single top production, whose 
cross section is proportional to |Vi;,p plus terms involving anomalous couplings. Hence, the complete 
determination of the Wtb coupling requires the combination of measurements in single top and top 
pair production. 

^The SM expectation, r{t Wb) ^ 1.6 GeV at the tree level, is one order of magnitude smaller 
than the width of the top invariant mass distribution reconstructed in the detector, which is about 
12 GeV (see for example Ref. [18]). Thus, deviations from the SM prediction for the top quark width 
are not likely to be observable. 
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3. The angular distribution of the W in the top rest frame, with respect to the top 
spin direction. 

The second class of observables, those related to W helicity fractions, may be defined 
(and in principle measured) for the decay of a top quark independently of the pro- 
duction mechanism, centre of mass energy, etc. In particular, these observables can 
be measured in single top as well as in top pair production. Previous literature [19] 
has already studied the sensitivity of the ATLAS experiment [20] for the measurement 
of helicity fractions in top pair production and subsequent semileptonic or dileptonic 
decay. Here we extend that analysis by including additional observables defined in 
Ref. [21]: the ratios of helicity fractions (denoted as "helicity ratios") and some new 
angular symmetries defined in the W rest frame. The dependence on anomalous cou- 
plings of the observables in the second class has been obtained in Ref. [21], including 
quadratic terms and keeping the b quark mass nonzero. In particular, it has been found 
that some of the new observables in Ref. [21] have a stronger parametric dependence 
on anomalous couplings than helicity fractions. In this work we study in detail the 
ATLAS sensitivity for their measurement in the semileptonic channel, paying a special 
attention to systematic uncertainties, both the theoretical ones and those related to 
the measurement in a real detector. We will eventually find that some of the new 
observables have smaller systematic uncertainties, and that their measurement can 
significantly improve the precision in the determination of anomalous couplings. 

This paper is organised as follows. In section [2] we briefiy set our notation and 
define the observables studied. A more extensive discussion of the theoretical aspects 
such as the relations among observables, their analytical expressions and plots of their 
dependence on anomalous couplings can be found in Ref. [21]. In section [3] the gener- 
ation of the ti signal and backgrounds is outlined, together with the selection criteria 
used to analyse them. In section [1] we present our results for the expected experi- 
mental measurement of the observables considered, and in section [5] we discuss their 
implications for the experimental determination of Wtb anomalous couplings. Section 
M is devoted to our conclusions. 
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2 The effective Wtb vertex and angular distributions 
in W rest frame 



The most general Wtb vertex containing terms up to dimension five can be written as 

= -^bnVLPL + VitPR)tW- 

-jf ^ i9LPL + 9rPr) t W- + h.c. , (1) 

with q = Pt — Pb the W boson momentum. The new anomalous couplings Vr, Ql and 
qr [11,21] can be related to //^, and in Ref. [19] (and references therein) as 
fi — Vr; — ~9l and — —Qr- If we assume CP is conserved these couplings can 
be taken to be real. Within the SM, Vl = Vth — ^ and the other couplings vanish at 
the tree level, while small nonzero values are generated at one loop level in the SM [22] 
and its extensions (see for example Refs. [16, 17]). 

The measurement of angular distributions and asymmetries in top decays can only 
determine ratios of couplings. (Besides, a moderate deviation from Vi, ~ 1 is not 
visible in top pair production and decay, as long as the top quark mainly decays to 
W^h and all other channels are rare.) Then, the value of Vl sets the global scale for 
the measurement of Vr, ql and qr in top decays. In this work we will normalise Vl 
to unity, and the limits on anomalous couplings presented correspond to Vl = 1- For 
any other value, the corresponding limits on anomalous couplings can be obtained by 
multiplying by the new Vl- 

It must be noted that, apart from the direct measurement at LHC, low-energy 
measurements already set indirect limits on non-standard Wtb couplings. The size 
of a Vr term is constrained by the measured rate of Br(6 57) = (3.3 ± 0.4) x 
10~^ [23]. A right-handed coupling \Vr\ > 0.04 would in principle give a too large 
contribution to this decay [24] which, however, might be (partially) cancelled with 
other new physics contributions. Hence, the bound IVr] < 0.04 is model dependent 
and does not substitute a direct measurement of this coupling. For ql the limits from 
b ^ are of the same order, while for qr they are much looser [25]. Besides, if one 
allows all anomalous couplings to be nonzero, direct and indirect limits turn out to be 
complementary, because they constrain different combinations of anomalous couplings. 

As we have already pointed out, the polarisation of the W bosons produced in the 
top decay is sensitive to non-standard Wtb couplings [26]. W bosons can be produced 
with positive, negative or zero helicity, with corresponding partial widths Vr, Vl, Vq 
which depend on V^, Vr, ql and qr. (General expressions for Vr, Vl, Vq in terms 
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of these couplings can be found in Ref. [21].) Their absolute measurement is rather 
difficult, so it is convenient to consider instead the helicity fractions Fj = Fj/F, with 
r = F/j + Fl + Fo the total width for t Wb. Within the SM, Fq = 0.703, Fl = 0.297, 
Fr = 3.6 X 10"^ at the tree level, for rut = 175 GeV, Mw = 80.39 GeV, = 4.8 GeV. 
We note that Fr vanishes in the = limit because the b quarks produced in top 
decays have left-handed chirality, and for vanishing nib the helicity and chirality states 
coincide. These helicity fractions can be measured in leptonic decays W — ^ iu. Let 
us denote by 9} the angle between the charged lepton three-momentum in the W rest 
frame and the W momentum in the t rest frame. The normalised angular distribution 
of the charged lepton can be written as 

^ l{l + cose;YFn + l{l-cosd;fFL + ^sm'e;F,, (2) 



Fdcos^; 8' - ' gv - ■ 4 

with the three terms corresponding to the three helicity states and vanishing interfer- 
ence [27]. A fit to the cos 6*1 distribution allows to extract from experiment the values 
of Fj, which are not independent but satisfy Fr + F^ + Fq = 1. From these measure- 
ments one can constrain the anomalous couplings in Eq. ([T]). Alternatively, from this 
distribution one can measure the helicity ratios [21] 

Pr,l = -jr- - -^r- ' l-^j 

J- -TO 

which are independent quantities and take the tree-level values pr = 5.1 x 10^"^, pl = 
0.423 in the SM. As for the helicity fractions, the measurement of helicity ratios sets 
bounds on Vr, ql and qr. 

A third and simpler method to extract information about the Wtb vertex is through 
angular asymmetries involving the angle 9}. For any fixed z in the interval [—1, 1], one 
can define an asymmetry 

_ A^(cos 9} > z)- N{cos 9} < z) 
" ~ N{cos9* > z)+ N{cos9* < z) ' ^' 

The most obvious choice is z = 0, giving the forward-backward (FB) asymmetry A^-q 
[11,28]§| The FB asymmetry is related to the W helicity fractions by 

A^^ = -^Fr-Fj}^. (5) 

Other convenient choices are z = =f(2^/'^ — 1). Defining /? = 2^/^ — 1, we have 
^ = _(22/3_i) ^ A, = A+ = 3P[Fo + {1 + P)Fr\, 



z 



(22/3 _i) ^ A, = A. = -3(3[Fo + il + /3)FL]. (6) 



^Notice the difference in sign with respect to the definitions in Refs. [11,28], where the angle 
, = TT — 9} between the charged lepton and b quark is used. 
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Thus, (^-) only depend on Fq and Fr (-Fl). The SM tree-level values of these 
asymmetries are Ap^ = —0.2225, = 0.5482, A_ = —0.8397. They are very sensitive 
to anomalous Wtb interactions, and their measurement allows us to probe this vertex 
without the need of a fit to the cos6'| distribution. We also point out that with 
a measurement of two of these asymmetries the helicity fractions and ratios can be 
reconstructed. For instance, using Eqs. ([6]) and requiring Fr + Fl + Fq = 1, it is found 
that 

J^R — Z n + 



1-/5 3(3(1 -(3^) 
1 A+ - (3 A. 



'2\ 



1 - (3 3(3(1 - (3' 



3 Simulation of signals and backgrounds and event 
selection 

The tt W^bW~b events in which one of the W bosons decays hadronically and the 
other one in the leptonic channel W — > iui (with i = e^,/i^), are considered as signal 
events. (From now on, the W boson decaying hadronically and its parent top quark will 
be named as "hadronic", and the W decaying leptonically and its parent top quark will 
be called "leptonic".) Any other decay channel of the tt pair constitutes a background to 
this signal. Top pair production, as well as the background from single top production, 
is generated with TopReX 4.10 [29] with default settings. Further backgrounds without 
top quarks in the final state, i.e. bb, W+jets, Z/7*+jets, WW, ZZ and ZW production 
processes, are generated using PYTHIA 6.206 [30]. In all cases we use CTEQ5L parton 
distribution functions (PDFs) [31]. Events are hadronised using PYTHIA, taking also 
into account initial state radiation (ISR), final state radiation (FSR) and pile-up. 

The generated background and signal events are passed through the ATLAS fast 
simulation packages ATLFAST 2.53 [32] and ATLFASTB [32]. These packages simulate 
the energy deposition in the calorimeter cells of all the stable particles in each event. 
The calorimeter cells are clustered within a cone of AR = ^(A0)2 + (Ar/)2 = 0.4, 
with the azimuthal angle and t] the pseudorapidity. Cells with transverse energy 
Ft > 1.5 GeV are used as cluster seeds and the cone algorithm is applied in decreasing 
order of Ft- Only clusters with Ft > 5 GeV are considered. The polar angle and 
the momentum of photons are smeared according to Gaussian parameterisations. For 
electrons, their momenta are smeared according to a Gaussian parameterisations. The 
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momentum of each muon is smeared according to a resolution which depends on the 
transverse momentum Pt, as well as on \ri\ and 0. The photon (electron) energy 
resolution is 5E/E < 2.9% (3.3%), for E > 20 GeV. The transverse momentum 
resolution of muons with pt < 100 GeV is Spr/pr ^ 2%. Photons, electrons and 
muons are selected only if they have \ri\ < 2.5 and pr > 5 GeV {pr > 6 GeV for muons). 
They are classified as isolated if the transverse energy of the cluster associated to the 
particle, inside a cone of AR = 0.2, does not exceed the particle energy by 10 GeV, 
and the AR from other energy clusters must be above 0.4. The clusters of energy 
depositions not associated to isolated photons, electrons or muons are used for the 
jet reconstruction. Their momenta are smeared according to a Gaussian distribution 
which depends on \r]\. Jets are selected if they have Et > 10 GeV. For E > 20 GeV, 
the jet energy resolution is better than 12% (for pseudorapidities |7;| < 3) and better 
than 24% (for \ri\ > 3). The missing transverse momentum is estimated by summing 
the transverse momentum of the isolated photons, electrons, muons and jets. The 
non-isolated muons and the clusters of energy deposition which are not associated to 
isolated photons, electrons, muons or jets, are also taken into account. In the ATLAS 
detector, it will be possible to identify 6 jets with \r]\ < 2.5 by using b tagging tools. The 
algorithm was simulated by setting a 6-tagging efficiency to 60%, with contamination 
factors set to 14.9% and 1.1% for c jets and light jets, respectively (the latter from 
light quark, gluon and tan leptons). In order to check the dependence of the analysis 
with the 6-tagging efficiencies, different values, 50% and 70% (corresponding to the 
expected 6-tag variation within the interesting signal transverse momentum range), 
were also considered for the systematic studies, with contamination factors of 9.2% 
(0.4%) and 23.3% (2.9%) for c-jets (light jets). 

Due to the hadronisation and FSR, the jets are reconstructed with less energies 
than those from the original quarks or gluons. The jets energies are calibrated by 
the ATLFASTB package, by applying a calibration factor, XJ^* = which is 

the ratio between the true parton energy and the reconstructed jet energy, obtained 
from reference samples [32]. The calibration factor depends on pt and is different for 
6-tagged and light jets. 

Signal events have a final state topology characterised by one isolated lepton (the 
isolation criterium requires the absence of additional tracks with px > 10 GeV inside 
a cone of AR — 0.4 around the lepton direction), at least four jets (among which 
exactly two must be tagged as b jets) and large transverse missing energy. We apply 
a two-level probabilistic analysis, based on the construction of a discriminant variable 
which uses the full information of some kinematical properties of the event. In the 
first level (called the pre-selection), a cleaner sample is obtained accepting events with: 
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(i) exactly one charged lepton with > 25 GeV, \ri\ < 2.5; (ii) at least 4 jets with 
Pt > 20 GeV, \ri\ < 2.5, two of them tagged as 6 jets and at least two not 6-tagged; (iii) 
missing transverse momentum above 20 GeV. The number of signal and background 
events (normalised to L = 10 fb~^) and the signal efficiency after the pre-selection are 
shown in the first column of Table [H Distributions of relevant variables are presented 
in Fig. [U 



Process 


Pre-Selection 


Final Selection 


ti iubbqq' 


262111 (11%) 


220024 (9%) 


tt (other) 


36745 


27060 


Single t 


12410 


7600 


Z+jets 


566 


253 


IV+jets 


3627 


1307 


WW, zz, zw 


109 


51 


total SM bkg. 


53457 


36271 



Table 1: Number of signal tt — ^ iubbqq' and background events, normalised to L = 
10 fb~^, after the pre-selection and final selection. The bb background is negligible after 
selection. 

The hadronic W reconstruction is done from the two non-6 jets with highest trans- 
verse momentum. The invariant mass of these two jets is represented in Fig. [H at the 
pre-selection. The mass of the hadronic top , also shown in Fig. [U is reconstructed as 
the invariant mass of the hadronic W and the b jet (among the two with highest px) 
closer to the W . The leptonic W momentum cannot be directly reconstructed due to 
the presence of an undetected neutrino in the final state. Nevertheless, the neutrino 
four-momentum can be estimated by assuming the transverse missing energy to be the 
transverse neutrino momentum. Its longitudinal component can then be determined, 
with a quadratic ambiguity, by constraining the leptonic W mass (calculated as the 
invariant mass of the neutrino and the charged lepton) to its known on-shell value 
M^r ~ 80.4 GeV. In order to solve the twofold quadratic ambiguity in the longitudinal 
component it is required that the hadronic and the leptonic top quarks have the min- 
imum mass difference. The reconstructed mass of the leptonic top is shown in Fig. [T] 
at the pre-selection. 

In the second level (the final selection), for each event we construct signal and 
background-like probabilities, p^s"^^' and V^^^^' , respectively, using probability density 
functions (p.d.f.) built from relevant physical variables: 
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(f) 
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(g) 




400 AT>0 500 



(h) 




■",(iep,(GeV/c^) 



(i) 



Figure 1: Kinematical distributions at the pre-selection level for the transverse mo- 
mentum of the charged lepton (a), the neutrino (b), the pr of the two non b jets used 
in the hadronic W reconstruction (c),(d), the 6 jet from the hadronic (e) and leptonic 
(f) top quarks. Invariant mass distributions of the hadronic W boson (g), the hadronic 
top (h) and the leptonic top (i). The ti signal (full line) and the SM backgrounds 
(shaded region) are normalised to L = 10 fb~^. 
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• The hadronic W mass. 



• The hadronic and leptonic top masses. 

• The transverse momentum of the 6-jets associated to the hadronic and the lep- 
tonic top quarks. 

• The transverse momentum of the jets used in the hadronic W reconstruction. 

These seven variables are shown in Fig. [U (c-i). Signal (£5 = U'^^iV-^^"'"'^) and back- 
ground {Cb = I^^^iV^"'^^') likelihoods (with n = 7, the number of p.d.f.) are used 
to define a discriminant variable = \ogiQ Cs/'^b- This variable is shown in Fig. [2] 
for the signal and background. The final event selection is done by applying a cut 
Lr > —0.2 on the discriminant variable, which corresponds to the highest S/ VB ratio. 
The number of background events (normalised to L = 10 fb~^) and signal efficiency 
after the final selection are shown in the second column of Table [TJ 



Figure 2: Discriminant variable for the SM background (shaded region) and the tt 
signal (full line), normalised to L = 10 fb~^. 

Two more cut-based analyses, omitted here for brevity, have also been performed. 
The results obtained depend more on the top mass reconstruction method than on the 
type (cut-based or probabilistic) of analysis performed. A detailed comparison of the 
three of them can be found in Ref. [33], where it is shown that the probabilistic analysis 
presented here gives the best results, with smaller systematic uncertainties. 



45000 



X 4C0C0 
<L> 
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4 Experimental measurement of angular distributions 
and asymmetries 



The experimentally observed cos6'£ distribution, which includes the ti signal as well as 
the SM backgrounds, is affected by detector resolution, ti reconstruction and selection 
criteria. In order to recover the theoretical distribution, it is necessary to: (i) subtract 
the background; (ii) correct for the effects of the detector, reconstruction, etc. For 
this purpose, we use two different sets of signal and background event samples: one 
"experimental" set, which simulates a possible experimental result, and one "reference" 
set, which is used to parameterise the effects mentioned and correct the previous sam- 
ple. The procedure is as follows. After subtracting reference background samples, the 
"experimental" distribution is multiplied by a correction function fc in order to recover 
the theoretical one expected in the SmQ The correction function is calculated, for each 
bin of the cos6'£ distribution, dividing the number of events at the generator level by 
the number of events after the event selection, using the reference sample. The "exper- 
imental" cos 6*1 distribution obtained after the simulation is shown in Fig. [3l together 
with the correction function obtained from the reference sample. The asymmetries are 
measured with a simple counting of the number of events below and above a specific 
value of cos^^ as in Eq. ((61). The procedure to correct for detector and reconstruction 
effects is basically the same, but with the cos 6*^ distribution divided into two or three 
bins. This has the advantage that the asymmetry measurements are not biased by the 
extreme values of the angular distributions, where correction functions largely deviate 
from unity and special care is required (see Fig. [3]). The helicity fractions and ratios ob- 
tained from a fit to the corrected distribution, as well as the angular asymmetries Afb, 
A±, are collected in Tabled with their statistical uncertainties. For easy comparison, 
we also include the theoretical values obtained at the generator level. 

Due to the excellent statistics achievable at LHC, which is apparent in Table [2], 
systematic errors play a crucial role in the measurement of angular distributions and 
asymmetries for a luminosity of 10 fb~^ or larger. A thorough discussion of the different 
systematic uncertainties in the determination of the correction functions is therefore 
compulsory. We estimate the systematic errors in the observables studied (asymme- 
tries, helicity fractions and ratios) by calculating them with various reference samples 



"Correction functions are determined assuming that the charged lepton distribution corresponds 
to the SM one. In case that a deviation from SM predictions (corresponding to anomalous couphngs) 
is found, the correction function must be modified accordingly, and the theoretical distribution re- 
calculated in an iterative process. These issues have been analysed in detail in Ref. [19], where it is 
shown that this process quickly converges. 
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2bC0C 




(a) (b) 

Figure 3: Simulated cos^^ distribution (a) and its correction function (b). In the first 
plot the ti signal (full line) and the SM backgrounds (shaded region) are normalised 
to L = 10 fb"\ 







Fl 


Fr 


PL 


PR 








Th. 


0.703 


0.297 


3.6 X 10-^ 


0.423 


5.1 X 10^^ 


-0.2220 


0.5493 


-0.8402 


Rec. 


0.700 


0.299 


0.0006 


0.4274 


0.0004 


-0.2231 


0.5472 


-0.8387 


Astat. 


0.003 


0.003 


0.0012 


0.0080 


0.0021 


0.0035 


0.0032 


0.0018 



Table 2: Theoretical and reconstructed values of helicity fractions, helicity ratios and 
angular asymmetries, with their statistical errors for L = 10 fb~^. 



and observing the differences obtained. In some cases the estimates are conservative, 
and they are taken as a reference for better comparison with previous analyses [19]. 
We consider uncertainties originating from: 

Monte Carlo generator: The correction functions obtained from a sample generated 
with TopReX are applied to a sample generated with ALPGEN [34]. The difference 
between the values obtained at the generator level and after the simulation is considered 
as systematic uncertainty. 

Structure functions: The correction functions obtained from a reference sample 
generated with CTEQ5L PDFs are applied to samples generated with CTEQ6L and 
MRST2001 PDFs in order to estimate the effects on the correction functions, and 
thus on the observables. The most significant deviations found are considered as the 
systematic error associated to the structure functions. 

Top mass dependence: Samples corresponding to top masses of 170, 175 and 180 GeV 
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are simulated, and the influence of rrit on the values obtained for the corrected observ- 
ables (using correction functions for = 175 GeV) is determined. The systematic 
error used here is obtained from a linear fit of the values found corresponding to a top 
mass uncertainty of 2 GeV. 

ISR and FSR: Their effect is studied following Ref. [35]. An event sample is created 
in which ISR and FSR are switched off in the event simulation. We compare the 
results of the reference sample (with ISR and FSR) with those obtained adding to it 
a normalised fraction of the sample without ISR nor FSR (from 0% to 25%, in steps 
of 5%). The values obtained for the observables are fitted with a linear function and 
the systematic error is considered as the effect of the presence of 20% (a conservative 
estimate of our level of knowledge of ISR and FSR) of the sample without ISR, FSR. 

b jet tag efficiency: The value of the 6 jet tag efficiency (and the corresponding c jet 
and light jet rejection factors) is varied from 50 to 70%, in steps of 5%, and the values 
obtained for the observables are fitted with a linear function The systematic error is 
considered as the effect on the observables of a variation of 5% in the b jet tagging 
efficiency, as compared with the standard value of 60%. 

b jet energy scale: The value of the b jet energy scale is changed from -5 to +5%, 
and the values obtained for the observables are fitted with a linear function. The 
systematic error is considered as the effect of a variation of 3% in the b jet energy scale. 

Light jet energy scale: The value of the energy scale of the light jets is changed 
from -3 to +3%, and the values obtained for the observables are fitted with a linear 
function. The systematic error is considered as the effect of a variation of 1% in the 
energy scale of the light jets. 

Background: The background (as obtained from the reference sample) subtracted 
to the selected sample is varied from -25 to 25%, in steps of 5%, and the values obtained 
for the observables are fitted with a linear function. The systematic error is considered 
as the effect of a variation of 10% on the background level (which takes into account 
the uncertainties in the cross-sections). 

Pile-up: The effect of pile-up events (2.3 events in average) is studied by comparing 
the values of the observables obtained with and without adding pile-up events. 

b quark fragmentation: The parameter e?, in the Peterson parameterisation for b 
quark fragmentation is changed from -0.006 to -0.0035, and the values obtained for the 
observables compared. The difference is considered as systematic error [35]. 

The systematic errors in each observable, resulting from these theoretical and sim- 
ulation uncertainties, are collected in Table [3l It can be observed that pn and A_ have 
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very small total systematic errors. In the case of pr, the improvement over is due 
to the cancellation of some of the systematic errors in the ratio, while the opposite 
happens in the case of p^, compared to Fl. 



Source 




Fl 


Fr 


Pt 




^FB 




A_ 


MC generator 


0.0002 


0.0002 


0.0004 


0.0006 


0.0000 


0.0035 


0.0015 


0.0006 


PDFs 


0.0032 


0.0022 


0.0009 


0.0046 


0.0008 


0.0021 


0.0005 


0.0014 


Top mass 


0.0065 


0.0060 


0.0006 


0.0124 


0.0007 


0.0034 


0.0039 


0.0005 


ISR+FSR 


0.0116 


0.0113 


0.0003 


0.0218 


0.0001 


0.0046 


0.0049 


0.0011 


h tag eff. 

o 


0.0065 


0.0062 


0.0003 


0.0126 


0.0003 


0.0039 


0.0046 


0.0004 


Eb scale 


0.0028 


0.0030 


0.0002 


0.0061 


0.0002 


0.0021 


0.0017 


0.0005 


Ej scale 


0.0034 


0.0037 


0.0002 


0.0074 


0.0002 


0.0038 


0.0023 


0.0014 


Back. 


0.0001 


0.0000 


0.0000 


0.0001 


0.0000 


0.0001 


0.0000 


0.0001 


Pile-up 


0.0091 


0.0086 


0.0005 


0.0175 


0.0002 


0.0080 


0.0051 


0.0006 


h frag. 


0.0049 


0.0037 


0.0012 


0.0078 


0.0011 


0.0045 


0.0000 


0.0012 


Total Asys. 


0.0189 


0.0178 


0.0018 


0.0356 


0.0016 


0.0130 


0.0099 


0.0028 



Table 3: Sources of systematic errors in the determination of helicity fractions, helicity 
ratios and angular asymmetries. 



The reduction of systematic errors compared to previous analyses deserves an ex- 
planation. In this analysis the W helicity fractions and ratios are obtained by fitting 
the angular distribution from —0.99 to +0.99, and it should be stressed that a depen- 
dence of the systematic error with the range of the fit has been observed. If the fit 
is performed between —0.89 and +0.89, the systematic errors on Fq, F^ and are 
respectively 0.0206, 0.0188 and 0.0033 (in good agreement with the results of Ref. [19]). 
However, if the fit is performed in the range [—0.89, +0.99] the results are respectively 
0.0190, 0.0182 and 0.0017, still in good agreement with the values on Table [3 This 
implies that the correct reconstruction of the most extreme bins of the angular dis- 
tribution is of utmost importance in order to control the error associated to the W 
polarisation measurements, if the fitting method is used. In the case of the asymme- 
tries, for A± the smaller errors are due to the greater stability of these measurements, 
obtained by counting events, compared to observables obtained from a fit to the cos 6*1 
distribution. We point out that the selection of z for the definition of A± in Eqs. ([6]) 
has not been optimised in order to achieve smaller systematic errors. Instead, these 
asymmetries have been defined in a simple way which allows to reconstruct easily the 
helicity fractions, using Eqs. ([7|). The results of our simulation, including statistical 
and systematic uncertainties, are summarised in Table [H 
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Observable 




Result 




Fo 


0.700 


±0.003 (stat) 


±0.019 (sys) 


Fl 


0.299 


±0.003 (stat) 


±0.018 (sys) 


Fr 


0.0006 


±0.0012 (stat) 


±0.0018 (sys) 


PL 


0.4274 


±0.0080 (stat) 


±0.0356 (sys) 


PR 


0.0004 


±0.0021 (stat) 


±0.0016 (sys) 


Afb 


-0.2231 


±0.0035 (stat) 


±0.0130 (sys) 


A+ 


0.5472 


±0.0032 (stat) 


±0.0099 (sys) 




-0.8387 


±0.0018 (stat) 


±0.0028 (sys) 



Table 4: Summary of the results obtained from the simulation for the observables 
studied, including statistical and systematic uncertainties. 

5 Limits on anomalous couplings 

With the results obtained in the previous section, summarised in Table [H and the para- 
metric dependence of the observables on Vr, ql and qr implemented in the computer 
program TopFit [21], constraints on the latter can be set. Naively, to obtain the la 
limit on a coupling x = Vr, ql, Qr derived from the measurement of some observable 
O, one would simply find the values of x for which O deviates la from its central 
valued Nevertheless, due to the quadratic dependence of the observables on Vr and ql 
near the SM point Vr = gi = 0, this procedure leads to overcoverage of the obtained 
confidence intervals [21], because their p.d.f. is not Gaussian even if the p.d.f. of the 
observable O is. In order to obtain the limits on an anomalous coupling x, given by the 
measurement of an observable O, we determine the p.d.f. of x numerically, using the 
acceptance-rejection method: we iteratively (i) generate a random value (with uniform 
probability) Xi within a suitable interval; (ii) evaluate the probability of 0(xi), given 
by the p.d.f. of O; (iii) generate an independent random number rj (with uniform 
probability); and (iv) accept the value Xi if the probability of 0(xj) is larger than rj. 
The resulting set of values {xj} is distributed according to the p.d.f. of x given by the 
measurement of O. The determination of a central interval with a given CL 7 is done 
numerically, requiring: (a) that it contains a fraction 7 of the total number of values 
{xj}; (b) that is central, i.e. fractions (1 — 7)/2 of the values generated are on each 
side of the interval. 

^This is the procedure originally followed in our previous work [33], as well as in Ref. [19]: For 
an observable O and a coupling x, intersecting the plot of 0[x) with the two horizontal lines O = 
Oexp ± AO, which correspond to the lu variation of O, gives the pretended la interval on x. 
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For X = this method gives results very similar to the intersection method in 
Refs. [19,33], whereas for Vr and gi the confidence intervals found are 20% and 30% 
smaller, respectively. The la limits derived from the measurement of each observable 
are collected in Table El assuming only one nonzero coupling at a time. We notice 
the improvement in sensitivity brought by the new observables Pr^l and A±: the best 
limits on Vr and ql are obtained from the measurement of pR, improving the limits 
from Fr by a factor of 1.13, and the best limits on qr are provided by ^4+, improving 
the limits from Fl by a factor of 1.34. This is due to the smaller (systematic plus 
statistical) uncertainties of these new observables and their stronger dependence on 
anomalous couplings. 





Vr 


9l 


9r 




{ql = gR = o) 


{Vr = gR = 0) 


{Vr = gL = 0) 






[-0.133,0.102] [ 


-0.0315,0.0219] 


Fl 


[-0.196,0.186] 


[-0.167,0.136] [ 


-0.0293,0.0212] 


Fr 


[-0.0373,0.1070] 


[-0.0491,0.0169] 




PL 


[-0.254,0.206] 


[ 


-0.0275,0.0227] 


PR 


[-0.0282,0.0987] 


[-0.0455,0.0129] 




^FB 


[-0.118,0.148] 


[-0.0902,0.0585] [ 


-0.0268,0.0227] 




[-0.140,0.146] 


[-0.112,0.0819] [ 


-0.0213,0.0164] 




[-0.0664,0.120] 


[-0.0620,0.0299] [ 


-0.0166,0.0282] 



Table 5: Limits on anomalous couplings obtained by the measurement of the observ- 
ables in the left column, with the constraint that only one non-standard coupling is 
allowed to be nonzero at a time. Dashes are shown where there is no significant sensi- 
tivity. 

These limits can be further improved by combining the measurements of the four 
observables pr^l and A±, including their correlations. We point out that the corre- 
lations among A±, p^ii do depend (as they must) on the method followed to extract 
these observables from experimental data. In our analysis A± are obtained by a simple 
event counting above and below a specific value of 2; = cos 6^1, while pr^l are obtained 
from a fit to the cos 6*1 distribution, divided in 20 bins. The correlations among these 
observables are derived as follows. We use a set of hypothetical "experimental mea- 
surements", in which each element of the set is a binned cos 9} distribution, as it would 
be experimentally obtained after correcting for detector effects. For each "measure- 
ment", the number of events in each cos 9} bin is obtained randomly using a Gaussian 
distribution centered at the expected SM value. We then calculate the average on this 
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set, denoted by (■), of the ten independent products of observables {A^A_), 
(A+pr), etc., where A±, p^^ are extracted from the 008 6*1 distribution as indicated 
aboveO The resulting correlation matrix is shown in Table [6l The correlations among 
A± and pr^l obtained are not affected by systematic uncertainties, as long as these do 
not significantly distort the shape of the cosO} distribution with respect to the SM one. 

A+ A_ PL PR 

A+ 1 0.1587 -0.8222 -0.1232 

A_ 0.1587 1 -0.08583 0.5688 

PL -0.8222 -0.08583 1 0.3957 

PR -0.1232 0.5688 0.3957 1 



Table 6: Correlation matrix for A±, Pr,l- 



When the four observables A± and pr^l are combined the assumption that only 
one coupling is nonzero can be relaxed. However, if Vr and ql are simultaneously 
allowed to be arbitrary, the limits on them are very loose and correlated, because for 
fine-tuned values of these couplings their effects on helicity fractions cancel to a large 
extent. In this way, values 0(0.4) of and ql are possible yielding minimal deviations 
on the observables studied. Therefore, in our combined limits, which are presented in 
Table [TJ we require that either Vr or ql vanishes. Limits for both Vr, ql nonzero 
require additional observables beyond the ones directly related to W helicity fractions, 
and will be presented elsewhere. 







Vr 


9l 


9r 




Pr,l 


[-0.0195,0.0906] 


X 


X 


A±, 


PR,L 


X 


[-0.0409,0.00926] 


X 


A±, 


Pr,l 


X 


X 


[-0.0112,0.0174] 


A±^ 


Pr,l 


X 


[-0.0412,0.00944] 


[-0.0108,0.0175] 


A±, 


Pr,l 


[-0.0199,0.0903] 


X 


[-0.0126,0.0164] 



Table 7: Limits on anomalous couplings obtained from the combined measurement of 
A±, Pr,l- In each case, the couplings which are fixed to be zero are denoted by a cross. 



^ Since the four observables A±, pu L are obtained from the corrected cos 6*1 distribution, there is no 
need to know the full kinematics of the ti event in order to determine their statistical correlation. On 
the other hand, if we are interested in, for example, the correlation between one of these observables 
and a top-antitop spin asymmetry, the full tt kinematics is needed. In the latter case, systematic 
errors can possibly influence the determination of the correlations. 
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For completeness, and to compare with previous literature we also present the 2a 
limits on non-standard couplings when only one of them is nonzero, 

Vr (2a) [-0.0566, 0.128] {ql = gn = 0) , 

gL (2a) [-0.0579, 0.0258] (Vr = gR = 0), 

gR (2a) [-0.0260, 0.0312] {Vr = gL = 0) . (8) 

A significant improvement, by factors of 3.25, 3.1 and 1.4, respectively, is obtained with 
the present analysis with respect to the results presented in Ref [19], which include the 
dilepton channel as well. This improvement is mainly due to: 

(i) The better sensitivity of the observables used. In the case of gL and Vr the 
improvement is moderate, with limits about 1.13 times smaller. For gR the 
improvement is more significant, by a factor of 1.34. 

(ii) The combination of pr^l and A±. 

(iii) The different statistical analysis used. For Vr and gL, the Monte Carlo method 
used to obtain the true 68.3% CL intervals also reduces their size by 20%-30%, 
as explained above. 

Finally, with the same procedure we obtain the 68.3% CL confidence regions on the 
anomalous couplings, presented in Fig.Hl The boundary of the regions has been chosen 
as a contour of constant x^- Iii case that the p.d.f. of Vr and gL were Gaussian, the 
boundaries would be ellipses corresponding to = 2.30 (see for instance Ref. [36]). In 
our non-Gaussian case the for which the confidence regions have 68.3% probability 
is determined numerically, and it is approximately 1.83 for the {gL^gn) plot and 1.85 
for {VR,gR). 

6 Conclusions 

In this paper we have investigated the ATLAS sensitivity to non-standard Wtb cou- 
plings. We have considered several observables: the helicity fractions Fj, helicity ratios 
Pr^l and angular asymmetries ApB, Although these observables can be defined 
and measured for any top production process with decay t — > Wb iub, we have con- 
centrated on top pair production at LHC with semileptonic decay, with a large cross 
section and in which the reconstruction of the final state is relatively easy. 

Due to the excellent statistics available at LHC, the precision reached is determined 
by systematic uncertainties. We have performed a very detailed study of the latter, 
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0.02 



V- Pl.r 



lo limit on gj^ 
lo limit on g^^ 



(a) 



-0.15 -0.1 -0.05 0.05 0.1 0.15 



A±- Pl.r 



(b) 



la limit on Vj^ 
Ic limit on 



-0.2 -0.1 0.1 



Figure 4: 68.3% CL confidence regions on anomalous couplings: gi and gn, for Vr = 
(a); Vr and gn, for gL = (b). The la combined limits in Table[7]are also displayed. 



both theoretical ones and from the experimental reconstruction. It has been found that, 
although the observables considered are theoretically equivalent (as noted in section [2]), 
the systematic uncertainties in the measurement of some of them, namely pn and 
are smaller. Since these observables also depend more strongly on anomalous couplings, 
their measurement provides a more sensitive probe for anomalous Wtb couplings than 
helicity fractions. Moreover, when the four measurements of pn^i and A± are combined, 
the sensitivity is further enhanced, reaching the 5.5%, 2.5% and 1.4% level for Vr, gi 
and gR in Eq. ([T]), respectively. This is an important achievement for a hadronic 
machine. Combining this measurement in tt semileptonic decays with the dilepton 
decay channel tt — ^ i~^iybi'~iyb and single top production will improve (to what extent 
is yet to be determined) these limits. 

Although providing probably the strongest limits, the observables studied in this pa- 
per are not sufficient to fully constrain anomalous Wtb couplings in a model-independent 
way. For nonzero Vr and gi, even of order 0(0.4), there are fine-tuned combinations for 
which their effects on helicity fractions and related observables almost cancel. Setting 
simultaneous limits on them requires additional observables with a different functional 
dependence on the Wtb couplings. For example, in the dilepton channel two spin 
asymmetries An' and Aee' involving the two leptons are found to be sensitive to Vr 
but rather independent of gL [21]. Ref. [19] has shown that these asymmetries can be 
measured with a good precision, 7% and 5%, respectively, and their study seems very 
promising. Spin asymmetries involving b quarks like A^f, and A£h exhibit a stronger 
dependence on anomalous couplings and when the appropriate detailed simulations 
are in place they will be studied. In addition single top production, involving Wtb 
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interactions in the production and the decay of the top quark will be studied, since 
it can provide complementary information about non-standard couplings through the 
cross sections for the different final states tj, ij, tb and ib, and spin asymmetries. 
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